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ABSTRACT 

We report the result of the analysis of the light curve of the microlensing event MOA-2009-BLG-016. The 
light curve is characterized by a short-duration anomaly near the peak and an overall asymmetry. We find that 
the peak anomaly is due to a binary companion to the primary lens and the asymmetry of the light curve is 
explained by the parallax effect caused by the acceleration of the observer over the course of the event due to 
the orbital motion of the Earth around the Sun. In addition, we detect evidence for the effect of the finite size 
of the source near the peak of the event, which allows us to measure the angular Einstein radius of the lens 
system. The Einstein radius combined with the microlens parallax allows us to determine the total mass of 
the lens and the distance to the lens. We identify three distinct classes of degenerate solutions for the binary 
lens parameters, where two are manifestations of the previously identified degeneracies of close/wide binaries 
and positive/negative impact parameters, while the third class is caused by the symmetric cycloid shape of the 
caustic. We find that, for the best-fit solution, the estimated mass of the lower-mass component of the binary is 
(0.04 ± 0.01) Mq, implying a brown-dwarf companion. However, there exists a solution that is worse only by 
A% 2 ~ 3 for which the mass of the secondary is above the hydrogen-burning limit. Unfortunately, resolving 
these two degenerate solutions will be difficult as the relative lens-source proper motions for both are similar 
and small (~ 1 mas yr -1 ) and thus the lens will remain blended with the source for the next several decades. 
Subject headings: gravitational lensing 
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For general microlensing events, the only lensing parameter 
that provides information about the physical parameters of the 
lens is the Einstein time scale fa- However, the time scale 
results from the combination of the underlying physical lens 
parameters of the lens mass M, relative lens-source parallax 
7r re i = AU (D^ 1 — Dg 1 ), and proper motion fi by 



te = — ; f? E = ( K M 7 r 1 . el ) 1 / 2 , 



(1) 



where n = 4G/(c 2 AU), 9e represents the angular Einstein ra- 
dius, and Z\ and £>s are the distances to the lens and source 
star, respectively. As a result, it is difficult to uniquely deter- 
mine the physical lens parameters from the time scale alone. 

For complete determination of the physical parameters of 
a lens, it is required to measure both the lens parallax and 
angular Einstein radius. The microlens parallax is defined by 
the ratio of the Earth's orbit to the Einstein radius projected 
on the observer plane, Fe, i.e., 



AU 

7T E = — . 



(2) 



In general, parallaxes are measured from the slight deviation 
of the overall sh ape of the light cu rve from a symmetric stan- 
dard light curve ( Paczvris kil 1 9861) . which is caused by the de- 
viation of the relative lens-source motion from a rectilinear 
motion due to the orbital motion of the Earth around the Sun 
dGouldll 1992b . Parallaxes are usually measured for long time- 
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FIG. 1 . — Light curve of the microlensing event MOA-2009-BLG-016. Also presented are the best-fit curves with and without the parallax effect. Middle panel 
shows the residual from the best-fit model (both parallax and finite-source effect included). The bottom panel shows the residuals from the best-fit finite-source 
(red) and point-source (blue) models in the region around the perturbation. The portion of the light curve corresponding to the time span of the bottom panel is 
shown by a box in the upper panel. We note that the data points in the bottom panel are binned by 4 hours to better show the difference. 



scale events for which the parallax effect is maximized. Ein- 
stein radii, on the other hand, are generally measured from the 
deviation of the light curve caused by the finite size of source 
stars such as caustic-approaching events (lGouldl fl994). With 
the measured parallax and Einstein radius, the mass of the lens 
and the distance to the lens are uniquely determined by 

M= — , (3) 

K7T E 

and 



7r E 6> E + 7rs ' 

where its = AU/Ds is the parallax of the source star. 

Unfortunately, the conditions for the optimal measure- 
ments of the parallax and Einstein radius are different 
and thus the chance to completely determine the physical 
parameters of lenses by measuring both quantities is low. 
In the literature, we find a total 13 microlensing events for 
which the lens mas ses were dete rmined. Th ese include 
EROS -2000-BLG-5 dAn et al.l l2002l). sc26-2218 dSmith et alj 
I2003al) . OGLE- 2002-BLG-018lKubas et alJ 120051). OGLE- 
2003-BLG-235 dBennett e t al.l|2006l), OG LE-2003-BLG-238 
dJiang et alj2004l). O GLE-2006-BLG- 1 09 dGaudi et alj2008t 
Bennett et alT "l2"009). OGLE-2 007-BLG-050 dBatista et all 
20091), OGLE-2 007-BLG-192 dBennett et all |2008|) OGLE- 
2007-BLG - 224 dGould et alj |2009|), O GLE12008-BLG-279 
dYeeetal.1 |2QQ9h. MACHO-LMC-5 dAlcock et alJ I200TI: 
iGouldetalj 120041). OGL E-2003-BLG-175/MOA-2003- 
BLG-045 dGhosh et alj 120041) and OGLE-2005-BLG-07 1 
dUdalski et alJl2005l:lDong et al.ll2009l) . 

In this paper, we report the result of the analysis of the 
microlensing event MOA-2009-BLG-016. We determine the 



mass of the lens and distance to it by measuring both the Ein- 
stein radius and lens parallax. We identify three distinct types 
of degeneracy, where two are previously known and the other 
is newly identified in this work. 

2. OBSERVATION 

The microlensing event MOA-2009-BLG- 
016 occurred on a star located at (RA,DEC) = 
(17°57'32.08",-34 h 21 m 10.06 s ), which corresponds to 
the Galactic coordinates of (l,b) = (-4.1°, -3.7°). It was first 
detected by the MOA collaboration on 2009 February 18 by 
using the 1.8 m telescope of Mt. John Observatory in New 
Zealand. An anomaly was detected on April 2 and it was 
announced to the microlensing community. In response to 
the alert, the /jFUN team conducted follow-up observations 
by using the 1.3 m SMARTS telescope of CTIO in Chile. 
The CTIO data are composed of 29 images in / band and 
5 images in V band. Data set from additional observatories 
were either single-epoch or near baseline and hence have not 
been included for analysis. Photometric reductions of the 
data were carried out by using the codes developed by the 
individual groups. 

In Figure Q] we present the light curve of the event. We 
note that the magnitude of the light curve is not calibrated. 
This is due to the lack of calibrated comparison stars in the 
field due to severe blending. However, the lensed star can 
be constrained from modeling combined with the color infor- 
mation obtained from the position of the source star on the 
color-magnitude diagram relative to the position of the center 
of clump giants in the field for which the de-reddened mag- 
nitude and color are well known. See more details in section 
3.1. We also note that data set from different observatories 
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and filters are aligned by fitting them a common model. 

3. CHARACTERIZATION 
3.1. Modeling 

Modeling microlensing light curves requires to include var- 
ious parameters. To describe light curves of standard single- 
lens events, a set of three parameters are needed. These in- 
clude the Einstein time scale, f E , the time of the closest lens- 
source approach, to, and the lens-source separation normal- 
ized by the Einstein radius at the time of the closest approach, 
mo- If light curves exhibit binary-induced anomalies, an addi- 
tional set of parameters is needed. These parameters include 
the mass ratio between the lens components, q, the projected 
binary separation in units of the Einstein radius, s, and the 
angle of the source trajectory with respect to the binary axis, 
a. For many cases of binary events, the normalized source ra- 
dius, = 9+/ 6b, is needed to describe the lensing magnifica- 
tion whenever the angular radius of the source star, 6+, plays 
an important role such as in the vicinity of a caustic cross- 
ing. For some long time-scale events with parallax-induced 
deviations, it is required to include the parallax parameters 
7TE,jv and 7Te £, which are the components of the lens-parallax 
vector 7Te projected on the sky in the north and east celestial 
coordinates, respectively. The direction of this vector is that 
of the lens-source relative motion in the frame of the Earth at 
the peak of the event. Similar to the deviation by the parallax 
effect, the source trajectory can also be affected by the orbital 
motion of the source if it is composed of binary stars ('xal- 
larap' effect). Under the assumption of a circular orbit and a 
very faint binary companion, the xallarap effect is parameter- 
ized by the orbital period P, inclination i, and phase angle tp 
of the orbit. 

Modeling light curves of microlensing events is a diffi- 
cult task. The large number of parameters makes brute-force 
searches of solutions difficult. In addition, a simple down- 
hill approach in the complicated \ 2 surface often results in 
wrong solutions of local minima. We, therefore, use a hybrid 
approach where grid searches are conducted over the space of 
the parameters of s, q, and a and the remaining parameters are 
searched by letting them vary so that they result in minimum 
X 2 at each grid of s, q, and a parameters. We use a Markov 
Chain Monte Carlo method for \ 2 minimization. Once the 
X 2 minima of the individual grid points are determined, the 
best-fit model is obtained by comparing the x 2 values of the 
individual grids. We investigate degeneracy of the solutions 
by probing local minima that appear in the space of the grid 
parameters. 

The light curve of the event MOA-2009-BLG-016 is char- 
acterized by two important features. One is the anomaly near 
the peak of the light curve and the other is the asymmetry 
of the overall light curve. From modeling, we find that the 
anomaly near the peak is well explained by the central pertur- 
bation caused by a close/wide binary. For the best-fit model, 
the determined values of the projected separation and mass 
ratio between the binary components are 

i = 0.21 ±0.01; 4 = 0.33 ±0.02. (5) 

We also find that the asymmetry of the light curve can be ex- 
plained by the parallax effect. The determined values of the 
parallax parameters for the best-fit parallax model are 

7T E , £ = 0.230 ±0.010; n EJf = 0.108 ±0.005. (6) 

18 Photometric data used for analysis are available upon requests. 
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FIG. 2. — Distribution of A\ 2 with respect to the xallarap parameters. The 
upper panel shows the distribution as a function the binary-source orbital pe- 
riod, P, and the lower panel shows the distribution with respect to the orbital 
phase, if), and inclination, i at a fixed orbital period of P = 1 yr. The position 
marked by "X" represents the position of the best-fit xallarap solution. 

The improvement of the fit with the addition of the parallax 
effect is A% 2 = 3641 . The parallax interpretation is consistent 
with the long time scale of the event, which is f E ~ 135 days 
for the best-fit model. 

It is known that a p arallax signal can b e mimicked by that 
of the xallarap effect (Smith et al. 2003b). To check this pos- 
sibility, we search for xallarap solutions under the assumption 
that the binary source is in a circular orbit. Figure|2]shows the 
distributions of A% 2 as a function of the xallarap parameters 
obtained from modeling. 

We find that the best fit xallarap solution yields an improve- 
ment of Ax 2 = 29 for 3 additional degrees of freedom. Al- 
though this improvement is formally highly significant and 
certainly implies that the 3 additional xallarap parameters are 
capable of "responding" to systematic deviations in the data 
that are not fully captured by the parallax modeling, we will 
now argue that the parallax solutions is preferred and the ad- 
ditional systematic deviations "detected" by the xallarap pa- 
rameters are most likely due to other effects. 

The first point is that the parallax signal (A% 2 = 3641) is 
more than 100 times stronger than the additional signal from 
xallarap, and such 1 % systematic effects are quite common in 
microlensing events. These may be due to real effects, such 
as binary orbital motion, third bodies in the system, or other 
physical effects, or may simply be due to systematics in the 
data, which (because of the fleeting nature of the events) are 
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FIG. 3. — Contours of A\ in the spaces of the binary parameters (s,q) [left panel], the parallax parameters (7rEJV,7rEfi) [middle panel], and the normalized 
source radius and the Einstein time scale [right panel] in the neighborhood of the best-fit solution (close I- model). 



taken under an extremely wide range of conditions. System- 
atic effects at this level can never be definitively tracked down 
because the A\ 2 is too small to adequately characterize the 
real effects (such as the mass and separation of a putative third 
body). 

Second, if we consider the entire two-dimensional space of 
possible binary orientations, only 1/6 lie within A\ 2 < 29 of 
the minimum, i.e., as close as the parallax solution (see Figure 
|2| . If we further account for the fact that all events must have 
parallax at some level, whereas only half of sources have a 
binary companion of any mass, with period P > lyr (see [2), 
then the prior probability of a xallarap solution is only 1/12. 
If there were known to be no other systematic effects (real or 
instrumental) affecting this event, then this 1/12 probability 
would have meager weight against a A\ 2 = 29. But since such 
low-level systematics are in fact common, the 1/12 probability 
must be taken seriously. In brief, while we cannot rule out the 
xallarap solution, we judge it to be relatively unlikely and so 
adopt the parallax solution. 

We find that a finite-source model is preferred over a point- 
source model with A% 2 = 12. The amount of A\ 2 is not big, 
but we note that it is statistically significant considering that 
the signal of the finite-source effect lasts only a short period 
of time during the source's approach close to the cusp of the 
caustic. 

In Figure Q] we present the best-fit model curve on the top 
of the light curve. Also presented in the top panel is the model 
curve without the parallax effect. Middle panel shows the 
residual from the best-fit model (both the parallax and finite- 
source effect included). The bottom panel shows the residu- 
als from the best-fit finite-source (red) and point-source (blue) 
models in the region around the perturbation. We note that the 
data points in the bottom panel are binned by 4 hours to better 
show the difference. 

In Figure [3] we present the contours of A\ 2 in the spaces 
of the binary parameters (s, q) [left panel], the parallax param- 
eters (7rE,£,7TE,Af) [middle panel], and the normalized source 
radius and the Einstein time scale (p*,?e) [right panel] in the 
neighborhood of the best-fit model. 

3.2. Degeneracy 

Although the basic characteristics of the lens system are de- 
fined, we find that there exist degenerate solutions. The upper 
panel of Figure [4] shows the local minima in the parameter 
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FIG. 4. — Likelihood contours for binary solutions as a function of the pro- 
jected binary separation (normalized by the Einstein radius) and the source 
trajectory angle (upper panel). Panels in the lower half show the geometry of 
the source trajectory (straight line with an arrow) with respect to the position 
of the caustic for the solutions corresponding to the individual local minima. 
The small yellow circle on the source trajectory represents the source star at 
the moment of the closest caustic approach. Among the 8 pairs of the local 
minima, we label only 4 pairs for which the values of A\ 2 from that of the 
best-fit model are relatively small. 

space of the projected binary separation and the source tra- 
jectory angle. Panels in the lower half show the geometry of 
the source trajectory with respect to the caustic for the solu- 
tions corresponding to the individual local minima. We note 
that among the total 8 pairs of local minima we label only 
4 pairs for which the values of A\ 2 from the best-fit model 
are relatively small. From the analysis of the individual local 
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TABLE 1 
Fit Parameters 



parameters 


close 1+ 


close I— 


wide 1+ 


wide I— 


close 11+ 


close II— 


wide 11+ 


wide II— 


X /dot 


649.1/645 


642.8/645 


655.9/645 


645.7/645 


663.7/645 


669.6/645 


673.0/645 


673.8/645 


s 


0.211 


0.208 


7.190 


7.296 


0.210 


0.207 


9.099 


8.934 


1 


0.298 


0.333 


0.686 


0.735 


0.582 


0.610 


5.157 


4.002 


a (deg) 


163.20 


-162.56 


163.53 


-162.91 


70.12 


-70.11 


69.22 


-69.18 


k (HJD') 


4923.805 


4923.779 


4923.914 


4923.882 


4923.928 


4923.934 


4923.780 


4923.785 


"0 


0.027 


-0.028 


0.021(0.027) 


-0.020(-0.027) 


0.036 


-0.036 


0.015(0.036) 


-0.016(-0.036) 


<E 


135.77 


135.29 


173.21(133.40) 


177.17(134.50) 


106.24 


108.55 


274.62(110.68) 


247.69(110.75) 


P* 


0.0055 


0.0056 


0.0046(0.0059) 


0.0043(0.0057) 


0.0087 


0.0085 


0.0031(0.0076) 


0.0035(0.0079) 


TE,N 


0.090 


0.108 


0.065(0.084) 


0.082(0.108) 


-0.026 


-0.023 


-0.011(-0.028) 


-0.010(-0.021) 


TTE.E 


0.224 


0.230 


0.171(0.221) 


0.180(0.237) 


0.275 


0.268 


0.104(0.257) 


0.117(0.261) 


fb,I 


0.692 


0.688 


0.690 


0.692 


0.594 


0.604 


0.600 


0.602 


fb.y 


0.693 


0.694 


0.690 


0.697 


0.604 


0.613 


0.604 


0.608 



NOTE. — HJD' = HJD -2450000. The parameters of the best-fit solution are marked in bold fonts. The parameters in the parentheses for 
wide solutions represent the values with respect to the mass of the binary component associated with the caustic involved with the perturbation. 
We note that fi, v and fi, ; represent the blended light fractions in V and / passbands, respectively, measured from the CTIO data. 



minima, we find that they result from three distinct types of 
degeneracy. 

The first ty pe of degenerac y is the well-known clo s e/wide 
degeneracy dDominik et al.1 1 1 999t lAlbrow et al.l 1 19991: 
lAfonso et all 120001: lAlbrow et al.1 120021k This degeneracy 
is caused by the similarity in shape between the caustics 
induced by a wide-separation binary with s > 1 and a close- 
separation binary with s < 1 . This can be seen in Figure [4] 
where one finds pairs of local minima with separations s > 1 
and s < 1 . 

The second type of degeneracy is caused by the mirror sym- 
metry between the pair of source trajectories with the im- 
pact para meters and source t rajectory angles of (uo,a) and 
(—uo,—a) (S mith e t al. 2003b). For a rectilinear motion, the 
two light curves resulting from the two source trajectories are 
identical. If the parallax effect is not negligible, however, the 
light curves from the trajectories are slightly different due to 
the curvature of the source trajectory. The pairs of local min- 
ima with trajectory angles a and -a in Figure|4]are caused by 
this degeneracy. 

The third type of degeneracy is caused by the shape of the 
caustic. When the caustic is produced by a binary with a sep- 
aration significantly larger or smaller than the Einstein radius, 
its shape is a symmetric cycloid with four cusps where two 
cusps are on the binary axis and the others are off the axis. 
Then, source trajectories approaching the caustic with angles 
a, a+ir/2, a+ir, and a+37r/2 result in a similar perturbation. 
We refer to this degeneracy as 'cycloid degeneracy'. We find 
that the caustic responsible for the central perturbation of the 
event MOA-2009-BLG-016 is very symmetric and thus the 
light curve is subject to this degeneracy. The four local min- 
ima on each quadrant of (s, a) parameter space in Figure|3]are 
caused by this degeneracy. 

In Table Q] we list the lensing parameters of the local min- 
ima along with values of \ 2 - We find that the models with 
trajectory angle of \a\ ~ 163° provide better fits than the cor- 
responding models with \a\ ~ 70°. Among the close-wide 
pairs of solutions with \a\ ~ 163°, we find close-binary mod- 
els are preferred. Among the two close-binary models with 
\a\ ~ 163°, we find that the "—Mo" model (close I-) provides 
the best fit to the observed light curve. 

3.3. Physical Parameters 

To determine we first determine the de-reddened mag- 
nitude /() and color (V - /)o of the source star from the mea- 
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FIG. 5. — Position of the source (lensed) star with respect to the centroid of 
clump giants in the instrumental (uncalibrated) color-magnitude diagram. 



sured offset between the source and the centroid of the clump 
giants in the instrumental color-magnitude diagram (CMD) 
constructed by using the V and I band images taken from 
CTIO (Figure|5J under the assumption that the source star and 
clump giants experience the same amount of extinction. The 
locations of the source star and the centroid of clump giants 
on the instrumental CMD are (V-I,I) S = (-0.86, 18.79) and 
(V-I,I) C = (-0.70, 16.70), respectively. The location of the 
source on the CMD is determined based on the light fractions 
of the source star in the / and V bands determined from mod- 
eling. With the known de-reddened magnitude and color of 
bulge clump giants of [(V-/)o,/o]c = (1-04, 14.32) toward the 
field, the de-reddened brightness and color of the source star 
are determined by (V-/)o.s = [(V-/) S -(V-7) C ] + 1.05 = 0.89 
and To.s = (h-Ic)+ 14.52 = 16.41, respectively. Here we 
adopt a mean dis tance to clump giants to ward the field of 8.8 
kpc estimated by Rattenb ury et al.l (120071) . Then, the angular 
source size is determined by first transformin g from (V - 7)p 
to (V - K)o using the color-color relation of iBessel & Brettl 
(1998) and then applying t he relation between (V - K)p and 
the angular stellar radius of Kervell a et al.l ( 120041) . The result- 
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TABLE 2 
Physical Parameters 



model 


Ax 2 


0E 








D L 


M 


Mi 


M 2 






(mas) 


(mas yr~' ) 


(mas yr~' ) 


(deg) 


(kpc) 


(Mq) 


(Mo) 


(M ) 


close 1+ 


6.3 


0.36±0.07 


0.97±0.17 


U2±Q.20 


65.9 


4.80±0.60 


0.18±0.03 


0.14±0.03 


0.04±0.()1 


close I- 


0.0 


0.36±0.06 


0.96±0.17 


1.13±0.20 


63.1 


4.70±0.58 


0.17±0.03 


0.13±0.02 


0.04±0.01 


wide 1+ 


13.1 


0.43±0.08 (O.33±0.O6) 


0.91±0.17 


1.05±0.20 


66.9 


4.98±0.64 


0.29±0.06 


0.17±0.03 


0.12±0.02 


wide I- 


2.9 


0.46±0.08 (O.35±0.06) 


0.94±0.16 


1.11±0.19 


63.6 


4.72±0.59 


0.28±0.05 


0.16±0.03 


0.12±0.02 


close 11+ 


21.0 


O.27±0.03 


0.91±0.10 


1.02±0.11 


90.1 


5.12±0.58 


0.12±0.01 


0.08±0.01 


0.04±0.01 


close II- 


26.8 


0.27±0.03 


0.90±0.11 


1.00±0.12 


89.7 


5.16±0.60 


0.12±0.02 


0.08±0.01 


0.05±0.01 


wide 11+ 


30.3 


0.74±0.14(0.30±0.06) 


0.99±0.19 


1.10±0.21 


90.9 


5.01±0.67 


0.88±0.19 


0.14±0.03 


0.73±0.16 


wide II- 


31.0 


0.64±0.11 (0.29±0.05) 


0.95±0.16 


1.06±0.18 


89.6 


5.07±0.65 


0.68±0.13 


0.14±0.03 


0.54±0.11 



NOTE. — Here 8e represents the angular Einstein radius, fj, = 0e/'e and fiQ represent the geocentric and heliocentric lens-source 
proper motion, respectively, ip^ is the angle of the proper motion with respect to the east, Dl is the distance to the lens, M is the total 
mass of the binary lens, and My and Mi are the masses of the binary components. We note that the subscript "1" is used to denote the lens 
component located closer to the source trajectory and thus My can be smaller than Mo . The parameters of the best-fit solution are marked 
in bold fonts. The values of A\~ are with respect to the best-fit solution, i.e. close I- model. The Einstein radius in parenthesis represents 
the value with respect to the mass of the binary component associated with the caustic involved with the perturbation. 



ing angular radius is 



% = (2.00 ±0.20) /xas, 



(7) 



where the uncertainty is estimated from the combination of 
the uncertainties of the colors and magnitudes of the source 
star and and an additional 7% intrinsic error in the c onversion 
proce ss from the measured color to source radius (lYee et al.l 
2009). With the measured source radius, the Einstein radius 
and lens-source proper motion are determined by 



and 



6» E (close I-) = 6J Pir = (0.36 ± 0.06) mas (8) 



^(close I-) = E /t E = (0.96 ±0.17) mas yr" 1 , (9) 



respectively. 

With the measured parallax and Einstein radius, the mass of 
the lens system and distance to the lens are determined from 
the relations in equation (3) and (4). For the best-fit model, 
these values are 



and 



M(close I-) = (0. 17 ± 0.03) M Q 



D L (close I-) = (4.70 ± 0.58) kpc, 



(10) 



(11) 



respectively. For the estimation of the uncertainty of Dl, we 
consider an 17% fractional error of the source location esti- 
mated by the bulge mass distribution model of lHan & Gould! 
(11995b . The values for other solutions are presented in Table 
|2] We note that, for the best solution, the estimated mass of 
the lower-mass component of the binary is 



M 2 (close I-) = 0.04 ± 0.01 M, 



0i 



(12) 



making it a brown-dwarf candidate. We note, however, that 
the estimated mass of the companion for the second-best 
model is above the hydrogen-buring limit. 

3.4. Resolution of Degeneracies 

It is found that analysis of the light curve of MOA-2009- 
BLG-016 alone results in degenerate solutions. However, it 



might be possible to resolve the degeneracies with extra in- 
formation. We check the possibility of resolving the degen- 
eracies from the measurement of proper motion from high- 
resolution observations. For this, we compute the heliocentric 
lens-source proper motion /j Q for the individual solutions. 

In Table Q] we present the magnitudes and directions of the 
proper motion vectors for the individual solutions. From the 
table, it is found that the high-resolution observation would 
be of limited use due to several reasons. First, the two sets of 
solutions caused by the cycloid degeneracy result in different 
directions of proper motion and thus resolution of the lens and 
source would make it possible to resolve the degeneracy, but 
this degeneracy is already clearly lifted with significant A% 2 
from the light curve alone. Moreover, this would be only pos- 
sible many years after the event considering the small magni- 
tude of the proper motion of P q ~ 1 mas yr" 1 . Second, the so- 
lutions caused by other degeneracies result in proper motions 
with not only a similar magnitude but also a similar direction, 
implying that the degeneracies would be difficult to be lifted 
even with high-resolution observations. 

4. CONCLUSION 

We analyzed the light curve of a microlensing event MOA- 
2009-BLG-016, which is characterized by high-magnification 
with an anomaly near the peak and an overall asymmetric 
light curve. We found that the anomaly and asymmetry of the 
light curve are explained by the lens binarity and the parallax 
effect, respectively. With the Einstein radius measured from 
the central perturbation combined with the lens parallax mea- 
sured from the overall asymmetric light curve, we determined 
the mass of the lens and distance to the lens. We identified 
three distinct types of degeneracy, where two were previously 
known and the other is first identified in this work. We also 
found that, for the best solution, the estimated mass of the 
lower-mass component of the binary is in the mass range of 
brown dwarfs. 
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